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I. INTRODUCTION AND OUTLINE OF THE MODEL 



During the last decade a large quantity of data on light nuclei with 4 < Z < 12 has 
become available. In particular, measurements of interaction cross sections by using ra- 
dioactive nuclear beams at intermediate and relativistic energies, have provided important 
data on nuclear radii The nuclear radius is a fundamental quantity which, in prin- 

ciple, provides information on the effective nuclear potential, shell effects and ground-state 
deformation. For exotic nuclei with extreme values of neutron to proton ratio, particularly 
important is the isospin dependence of nuclear radii which can signal the onset of new phe- 
nomena like, for example, the formation of skin and halo structures. Data on ground-state 
deformations are also very important for the study of shell effects in exotic nuclei. In partic- 
ular, they reflect major modifications in the shell structures, the disappearance of standard 
and the occurrence of new magic numbers. Different ground-state deformations of proton 
and neutron density distributions are expected in some nuclei with extreme isospin projec- 
tion quantum number T z . Another source of information on the effective nuclear potential 
in exotic systems at the limits of stability are the single-nucleon separation energies. The 
neutron drip-line has been reached for nuclei up to Z = 9 ||. On the proton-rich side the 
drip-line has been experimentally fully mapped up to Z=21, and possibly for odd-Z nuclei 
up to In || . In very neutron- rich nuclei the weak binding of the outermost neutrons causes 
the formation of the neutron skin on the surface of a nucleus, and the formation of one- 
and two-neutron halo structures. The established two-neutron halo nuclei are 6 He, n Li, and 
14 Be, and the one-neutron halo nuclei are 11 Be and 19 C. Recent data present evidence for 
a one-neutron halo in 22 N, 23 and 24 F. The formation of the neutron skin is well established 
in the neutron-rich Na isotopes , and the related phenomenon of low- lying pygmy isovec- 
tor dipole resonances has recently been observed in O isotopes [|10| . On the proton-rich side 
evidence has been reported for the existence of a proton skin in 20 Mg | llfl , and a beautiful 
example of exotic decay modes is provided by the two-proton emitter 18 Ne |12| . 



In the present work the Relativistic Hartree-Bogoliubov (RHB) model is applied in the 
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analysis of ground-state properties of Be, B, C, N, F, Ne and Na isotopes. Based on the 
relativistic mean-field theory and on the Hartree-Fock-Bogoliubov framework, the RHB 
model provides a unified description of mean-field and pairing correlations. It has been 
successfully applied in the description of nuclear structure phenomena in exotic nuclei far 
from the valley of (3- stability and of the physics of the drip lines. On the neutron-rich 
side RHB studies include: the halo phenomenon in light nuclei Hl3 |, properties of light 



nuclei near the neutron-drip [|14], the reduction of the spin-orbit potential in nuclei with 



extreme isospin values ||15|| , the deformation and shape coexistence phenomena that result 



from the suppression of the spherical N=28 shell gap in neutron-rich nuclei [16|, properties 



of neutron-rich Ni and Sn isotopes ||17|| . In proton- rich nuclei the RHB model has been used 



to map the drip line from Z = 31 to Z = 73, and the phenomenon of ground-state proton 



radioactivity has been studied [IS -20]. In a very recent study of the isovector channel of 



the RHB model [21], a very good agreement with experimental data has been obtained for 



ground-state properties of nuclei that belong to the A = 20 isobaric sequence. 



In the framework of the relativistic mean field (RMF) approximation |EJ nucleons are 
described as point particles that move independently in the mean fields which originate 
from the nucleon-nucleon interaction. The theory is fully Lorentz invariant. Conditions of 
causality and Lorentz invariance impose that the interaction is mediated by the exchange of 
point-like effective mesons, which couple to the nucleons at local vertices. The single-nucleon 
dynamics is described by the Dirac equation 



j-ia • V + (3(m + g a a) + g w u Q + g p T 3 p° 3 + e- 



LJVa°\ ft = e*lH. (1) 



a, u>, and p are the meson fields, and A denotes the electromagnetic potential. g a g u , and g p 
are the corresponding coupling constants for the mesons to the nucleon. The lowest order 
of the quantum field theory is the mean-field approximation: the meson field operators are 
replaced by their expectation values. The sources of the meson fields are defined by the 
nucleon densities and currents. The ground state of a nucleus is described by the stationary 
self-consistent solution of the coupled system of Dirac and Klein-Gordon equations. 



In addition to the self-consistent mean-field potential, pairing correlations have to be in- 
cluded in order to describe ground-state properties of open-shell nuclei. In the framework of 
the relativistic Hartree-Bogoliubov model, the ground state of a nucleus |$ > is represented 
by the product of independent single-quasiparticle states. These states are eigenvectors of 
the generalized single-nucleon Hamiltonian which contains two average potentials: the self- 
consistent mean-field T which encloses all the long range particle-hole (ph) correlations, and 
a pairing field A which sums up the particle-particle (pp) correlations. In the Hartree ap- 
proximation for the self-consistent mean field, the relativistic Hartree-Bogoliubov equations 
read 

h D -m-\ A \ (U k (r)\ (U k {r)\ 

-A* -h D + m + \) \V k {r)J \V h (r))' 

where hp is the single-nucleon Dirac Hamiltonian ([!]), and m is the nucleon mass. The 
chemical potential A has to be determined by the particle number subsidiary condition in 
order that the expectation value of the particle number operator in the ground state equals 
the number of nucleons. A is the pairing field. The column vectors denote the quasi-particle 
spinors and E k are the quasi-particle energies. 

The RHB equations are solved self-consistently, with potentials determined in the mean- 
field approximation from solutions of Klein-Gordon equations for the meson fields. The 
Dirac-Hartree-Bogoliubov equations and the equations for the meson fields are solved by 
expanding the nucleon spinors £4(r) and V^(r), and the meson fields in terms of the eigen- 
functions of a deformed axially symmetric oscillator potential. The calculations for the 
present analysis have been performed by an expansion in 14 oscillator shells for the fermion 
fields, and 20 shells for the boson fields. A simple blocking procedure is used in the calcu- 
lation of odd-proton and/or odd- neutron systems. The blocking calculations are performed 
without breaking the time-reversal symmetry. A detailed description of the Relativistic 
Hartree-Bogoliubov model for deformed nuclei can be found in Ref. [[19 . 
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II. GROUND-STATE PROPERTIES OF DEFORMED LIGHT NUCLEI 



In parallel with the experimental work of the last decade, many theoretical analyses 
have been performed of the structure of nuclei in the mass region 10 < A < 30. Both 
microscopic mean-field and shell-model approaches, as well as various microscopic cluster 
models, have been used to study properties of ground and excited states of isotopic and 
isobaric sequences, and to describe specific structure phenomena in exotic nuclei. It has 
been shown that theoretical models reproduce the global trends of nuclear sizes and binding 
energies. However, special assumptions have to be made, or even new models have to be 
designed, in order to describe more exotic phenomena like, for example, the location of the 
neutron drip line in Oxygen, or the ground-state deformation of 32 Mg. 

In Ref. |14| we reported spherical RHB calculations of neutron-rich isotopes of N, O, 
F, Ne, Na and Mg. By using several standard RMF effective interactions, we analyzed 
the location of the neutron drip-line, the reduction of the spin-orbit interaction, rms radii, 
changes in surface properties, and the formation of neutron skins and of neutron halos. It 
was shown that, even without taking into account the deformation of the mean field, the 
RHB model correctly describes the global trends of the observed ground-state properties. 
The exception is, of course, the location of the neutron drip line in Oxygen, which none of 
the RMF effective interactions reproduces. 

In the study of Ref. |ST[ we performed deformed RHB calculations of ground-state proper- 



ties of nuclei that belong to the A = 20 isobaric sequence. The NL3 effective interaction |23 



was used for the mean-field Lagrangian, and pairing correlations were described by the pair- 



ing part of the finite range Gogny interaction D1S |24 |. This particular combination of 
effective forces in the ph and pp channels has been used in most of our recent applications 
of the RHB theory. RHB results for binding energies, neutron and proton ground-state 
density distributions, quadrupole deformations, nuclear matter radii, and proton radii were 
compared with available experimental data. The very good agreement with the observed 
ground-state properties as function of the isospin projection T z , led to the conclusion that 
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the isovector channel of the NL3 interaction is correctly parameterized and that this effec- 



tive force can be used to describe not only medium-heavy and heavy nuclei fT7j-|20"|| , but also 
properties of relatively light nuclei far from /3-stability. 

In the present work we apply the RHB model, with the NL3+D1S effective interaction, in 
the analysis of ground state properties of Be, B, C, N, F, Ne and Na isotopic sequences. We 
perform deformed RHB calculations and compare radii, separation energies and quadrupole 
deformations with available experimental data and with the predictions of the finite range 
droplet model (FRDM) |32j . Since the RHB equations are solved in the configuration space 
of harmonic oscillator basis states, for nuclei at the drip lines we do not expect an accurate 
description of properties that crucially depend on the spatial extension of the wave functions 
of the outermost nucleons, especially on the neutron-rich side. Thus, we do not attempt to 
describe radii of halo nuclei. We also do not repeat the spherical calculations of O isotopes, 



which can be found in Ref. 14 



In Fig. [1] we display the proton, neutron and matter radii, ground-state quadrupole de- 
formations, and one-neutron separation energies of Beryllium isotopes, calculated with the 
NL3 + Gogny D1S effective interaction. The RHB values are compared with the experimen- 
tal radii |!]|| and separation energies [ ^Sfl . The calculated matter radii reproduce the trend 
of the experimental values, except for the halo nucleus n Be [ 26|,|27 |. For the proton radii, on 
the other hand, the theoretical values are considerably lower than the experimental ones [|TJ , 
especially for 14 Be, though it should be noted that the experimental proton radius has a 
very large uncertainty. The calculated neutron skin in 14 Be is very large: r n — r p = 0.71 
fm, and also the deformations of the proton and neutron densities in this nucleus are very 
different. 14 Be has a large prolate deformation /3 2 — 0.36 and the ratio of neutron to pro- 
ton quadrupole moments is Q n /Qp = 3.6. The one-neutron separation energies are also in 
excellent agreement with the empirical values In particular, 13 Be is predicted to be 

unbound by 180 keV. 

In Fig. the same comparison is shown for neutron-rich Boron isotopes. The calculated 
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matter radii are in excellent agreement with the recent experimental data || for 17 B and 
19 B, while they are larger than the older experimental values QXJ] for 14 B and 15 B. Unlike in 
the case of Be, the calculated proton radii for these two nuclei agree well with the empirical 
values, but the theoretical neutron radii are much larger. The RHB calculation also predicts 
14 B and 15 B to be spherical in the ground state, while the heavier Boron isotopes are strongly 



prolate deformed. The separation energies agree with the empirical values ||25|| , and we note 
that both 16 B and 18 B are predicted to be neutron unbound. 

The calculated quantities which characterize the ground-states of neutron-rich Carbon 
isotopes are displayed in Fig. || The proton radius of 14 C is compared with the experi- 
mental value from Ref. |p8|| , and the matter radii are shown in comparison with very recent 
experimental data 0. The trend of the experimental matter radii is clearly reproduced by 
the RHB calculation. Of course, for the one-neutron halo nucleus 19 C [^9|] the calculation in 
the harmonic oscillator configuration space cannot reproduce the anomalous matter radius. 
The RHB model predictions for the one-neutron separation energies are in agreement with 
the empirical values ||25|| , though the theoretical values are slightly larger for the even-N 
isotopes. In particular, from an analysis of the angular distribution of the 18 C +n center of 
50], the neutron separation energy is determined to be 530 ± 130 keV, and the RHB 



mass 

calculated value is 510 keV. 14 C and 22 C are spherical, and all other neutron-rich Carbon 
isotopes, except 17 C, have oblate ground states. The prolate minima, however, are found at 
relatively low excitation energies, and in Fig. |3] we also plot the energy differences between 
prolate and oblate minima. 

The RHB results for Nitrogen isotopes are shown in Fig. 01. The proton radii are compared 



with the experimental values from Ref. J28| . The calculated matter radii reproduce the global 
trend of the experimental data J7j, but not the sudden increase of the radii at N=15. In 
fact, in the recent measurement of the interaction cross sections for 14_23 N, 16_24 0, and 
i8-26p on Carbon targets at relativistic energies 0, a sharp increase of matter radii at 
N=15 was observed in all three isotopic chains (see also Fig. Rf). The deduced matter radii 
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for 22 N, 23 0, and 24 F are much larger than those of their respective neighbors with one 
neutron less, and the deduced nucleon density distributions show a long neutron tail for 
these nuclei, comparable to that in n Be. It was therefore concluded that these three nuclei 
exhibit a one-neutron halo structure. Since they are spherical, the halo structure should 
result from the odd neutron occupying the 2s\/2 orbital. The absence of the centrifugal 
barrier for an s-orbital facilitates the formation of the long tail of the wave function, i.e. 
the halo structure. However, in Ref. @ it was also noted that the one-neutron separation 
energies for 22 N, 23 0, and 24 F are larger than 1 MeV (1.22 ± 0.22 MeV, 2.74 ± 0.12 MeV 
and 3.86 ±0.11 MeV, respectively), and much larger than in n Be and 19 C. In a very recent 
analysis [^TJ it has been pointed out that the conventional fixed core-plus-neutron model for 
halo nuclei is unable to explain the observed enhanced cross section for these nuclei by any 
selection of neutron orbitals, and therefore a possibility of an enlarged core structure has 
been suggested. The mechanism of core enlargement in these three nuclei, however, has not 
been explained. We first note that the present RHB calculation reproduces the empirical 
one-neutron separation energies [ 25| . In particular, for 22 N we even obtain a slightly lower 
one-neutron separation energy, and the theoretical matter radius coincides with the one 
deduced from the experimental interaction cross section. Also for 23 N, the calculated and 
empirical separation energies coincide, and the theoretical matter radius is only slightly 
below the large experimental error bar. The main difference is in the matter radii of the 
lighter isotopes (a similar situation also occurs for the Fluorine isotopes, see Fig. §). The 
calculated radii are somewhat larger than the experimental values and therefore at N=15 do 
not display the sharp discontinuity which, in Ref. [0], is taken as evidence for the formation 
of the neutron halo. In the present calculation the gradual increase of matter radii reflects 
the formation of the neutron skin. This is shown in the upper right panel of Fig. f|, where 
the values of r n — r p are plotted as function of the neutron number, and in Fig. |^ we 
display the self-consistent RHB neutron and proton ground-state density distributions of 
even-N Nitrogen isotopes. It should be pointed out, however, that the formation of the halo 
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structure can only be observed if calculations were performed in coordinate space. Moreover, 
particle number projection might be necessary in order to reproduce the sharp increase of 
matter radii. Finally, we note that the RHB NL3+D1S calculation predicts the heaviest 
particle stable Nitrogen isotope to be 23 N, in excellent agreement with recent data on the 
neutron drip-line [§]. 

Very similar results are obtained for the Fluorine isotopes. In Fig. ^|we compare the RHB 
theoretical proton, neutron and matter radii, and one-neutron separation energies with the 
experimental radii an d separation energies [25]. The calculated matter radii do not 

reproduce the discontinuity at N=15, though for the heaviest isotopes they are found within 
the experimental error bars. We note, however, that in the case of Fluorine isotopes the 
calculated one-neutron separation energies do not reproduce the staggering between even-N 
and odd-N values below iV < 15. It is interesting to note that, like in the case of Nitrogen, 
the RHB model with the NL3+D1S effective interaction correctly predicts the location of 
the drip-line ||: the last bound isotope of Fluorine is 31 F. Therefore, in agreement with 
experimental data, we obtain that the neutron drip-line is at N=16 for Z=7, and at N=22 
for Z=9. On the other hand, none of the standard RMF effective interactions reproduces the 
location of the drip-line for Oxygen. It has been argued that the sudden change in stability 
from Oxygen to Fluorine may indicate the onset of deformation for the neutron-rich Fluorine 
isotopes ||. In the present calculation, however, all Fluorine isotopes up to 31 F turn out to 
be essentially spherical. 

In Ref. |13|| we performed spherical RHB calculations of the Ne isotopic chain. In par- 
ticular, we studied the formation of neutron halo structures in drip-line Ne nuclei (N > 20). 
It was shown that the properties of the lf-2p neutron orbitals near the Fermi level, and the 
neutron pairing interaction play a crucial role in the possible formation of the multi-neutron 
halo. In the present analysis we have performed calculations in the deformed harmonic os- 
cillator basis. In Fig. |7| we display the proton, neutron and matter radii, the ground-state 
quadrupole deformations and the one-neutron separation energies of Ne isotopes. The cal- 
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culated /?2 values are shown in comparison with the predictions of the finite-range droplet 
model |32)| , and the separation energies are compared with experimental data [25|| . The 
FRDM and the present RHB calculation predict a similar mass dependence of the ground- 
state quadrupole deformation. We note two spherical regions around A=16 and A=28. Both 
models reproduce the large prolate deformations around A=20, and predict prolate shapes 
in the region of possible halo structures A > 30. Pronounced differences in the predicted 
/?2 values are found for A=19 and A=24,25. The later probably indicates a region of shape 
coexistence. For A=28 the FRDM predicts an oblate (3 ~ —0.2 deformation, while a spheri- 
cal shape is calculated in the RHB model. The calculated separation energies reproduce the 
odd-even staggering and agree quite well with the experimental values. 

The ground-state properties of the Na isotopic sequence are illustrated in Fig. |[ The 
one-neutron separation energies are shown in comparison with experimental data [p5|1 . The 
calculated values reproduce the empirical staggering between even- and odd-A isotopes, 
although for A > 24 the theoretical separation energies are systematically somewhat larger 
for the even-N isotopes. The calculated radii are compared with the experimental data: 
matter radii ||, neutron radii ||, and charge isotope shifts pfl . An excellent agreement 
between theory and experiment is obtained. For the matter and neutron radii the only 
significant difference is at A=22, but this dip in the experimental sequence has recently 
been attributed to an admixture of the isomeric state in the beam ||. Except for the 
lightest isotope shown, i.e. 20 Na, the calculated charge isotope shifts reproduce the empirical 
A-dependence. A significant difference between the theoretical and experimental values is 
observed only for A > 29. The calculated ground-state quadrupole deformations of the 
Na isotopes are compared with the predictions of the finite-range droplet model |32[]. We 
note that, while the FRDM predicts all Na isotopes with A < 28 to be strongly prolate 
deformed, the result of the RHB calculation is the staggering between prolate and oblate 
shapes, indicating the onset of shape coexistence. In particular, 26 ' 27 Na are predicted to 
be oblate, while prolate ground-state deformations are calculated for 28,29 Na. Very recent 
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experimental data on quadrupole moments of 26_29 Na |34| confirm this prediction. 

Our results for matter radii are summarized in Fig. || where we plot the calculated values 
for the A = 14, 15, 16, 17, 18, 19 isobaric chains as functions of the isospin projection T z . An 
excellent overall agreement is found between the experimental data and the matter radii 
calculated with the NL3+D1S RHB effective interaction. 



Finally, in Fig. [10] the matter radii of mirror nuclei are compared as function of the 
isospin projection T z . It is interesting to note that for AT Z = 1, the nuclei with T z > have 
larger radii than their mirror partners with T z < 0. Due to the strong effect of the Coulomb 
interaction, on the other hand, for AT Z = 2 and 3 the proton-rich nuclei have almost 0.1 fm 
larger nuclei than their T z > mirror partners. 



III. SUMMARY 

This work presents an analysis of ground-state properties of Be, B, C, N, F, Ne and Na 
isotopes in the framework of the Relativistic Hartree-Bogoliubov (RHB) model. In the last 
couple of years this model has been very successfully applied in the description of nuclear 
structure phenomena in medium-heavy and heavy exotic nuclei far from the valley of (5- 
stability and of the physics of the drip lines. The present analysis covers a region which is 
probably at the low-mass limit of applicability of the mean-field framework. This work is 
also a continuation of our previous applications of the RHB model of Ref. |14| (spherical 
RHB calculations of neutron-rich isotopes of N, O, F, Ne, Na and Mg), and of Ref. pl| 
(deformed RHB calculations of ground-state properties of the A = 20 isobaric sequence). 

The present calculation has been performed in the configuration space of a deformed 
harmonic oscillator basis states. The NL3 effective interaction has been used for the mean- 
field Lagrangian, and pairing correlations have been described by the pairing part of the 
finite range Gogny interaction D1S. The calculated neutron separation energies, quadrupole 
deformations, nuclear matter radii, and differences in radii of proton and neutron distribu- 
tions have been compared with very recent experimental data. Here we summarize the main 
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results: 

• For the neutron drip-line nucleus 14 Be the RHB calculation predicts a large pro- 
late deformation @ 2 = 0.36 and the ratio of neutron to proton quadrupole moments 
Qn/Qp = 3.6. 

• For the neutron-rich Boron isotopes, the calculated matter radii reproduce the recent 
experimental data || for 17 B and 19 B. This is an important result, since 19 B (T z = 9/2) 
has one of the largest N/Z values known at present in the low- mass region of the nuclear 
chart. 

• Although the present calculation, performed in the deformed harmonic oscillator con- 
figuration space, cannot reproduce the anomalous matter radius of the one-neutron 
halo nucleus 19 C, the neutron separation energy 510 keV is in excellent agreement 
with the experimental value 530 ±130 keV obtained from an analysis of the angular 
distribution of the 18 C +n center of mass ||30|| . A large oblate deformation (3 2 ~ —0.4 
is predicted for 19 C. 

• The RHB model with the NL3+D1S effective interaction predicts the location of the 
neutron drip-line in Nitrogen and Fluorine in agreement with recent experimental 
findings ||: the heaviest particle stable isotopes are 23 N and 31 F. The calculation, 
however, does not reproduce the location of the neutron drip line in Oxygen. 

• The calculated matter radii of the neutron-rich Nitrogen and Fluorine isotopes are in 
agreement with very recent experimental data J7|. However, the calculation does not 
reproduce the sudden increase of the radii at N=15, which was taken as evidence for 
the formation of the neutron halo. 

• For the neutron-rich Neon isotopes the RHB model predicts pronounced prolate de- 
formations in the region of possible halo structures (A > 30). 
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For the Na isotopic sequence the calculated radii are in excellent agreement with 
experimental data on matter radii ||, neutron radii ||, and charge isotope shifts p3[ . 
The calculated ground-state quadrupole deformations confirm the recent experimental 
data on quadrupole moments of 26_29 Na fl34|. 
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FIG. 1. Proton, neutron and matter radii, ground-state quadrupole deformations and 
one-neutron separation energies of Beryllium isotopes, calculated with the NL3 + Gogny D1S 
effective interaction. The theoretical values are compared with the experimental radii |5Q and 



separation energies [25]. 



FIG. 2. Same as in Fig. |], but for Boron isotopes. 

FIG. 3. Proton, neutron and matter radii, ground-state quadrupole deformations, one-neutron 
separation energies and energy differences between prolate and oblate minima for Carbon isotopes. 
The theoretical values are compared with the experimental radii [^jTj and separation energies [25]. 

FIG. 4. The RHB theoretical proton, neutron and matter radii, and one-neutron separation 
energies of Nitrogen isotopes, compared with the experimental radii p8|,^j7t and separation ener- 



gies [25]. 



FIG. 5. Self-consistent RHB neutron and proton density profiles of Nitrogen isotopes. 

FIG. 6. The RHB theoretical proton, neutron and matter radii, and one-neutron separation 
energies of Fluorine isotopes, compared with the experimental radii [p8|,p|,|7[| and separation ener- 



gies [25]. 



FIG. 7. Proton, neutron and matter radii, ground-state quadrupole deformations and 
one-neutron separation energies of Neon isotopes. The calculated fli values are displayed in com- 
parison with the predictions of the finite-range droplet model [32|, and the separation energies are 
compared with experimental data 



FIG. 8. One-neutron separation energies, radii, charge isotope shifts and ground-state 
quadrupole deformations of Sodium isotopes. The RHB calculated values are compared with the 
experimental data: one-neutron separation energies p!|] , matter radii ||, neutron radii 0, charge 



isotope shifts p3fl, and with the @2 values calculated in the finite-range droplet model [32]. 
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FIG. 9. The nuclear matter radii of the A = 14, 15, 16, 17, 18, 19 isobaric chains as functions 
of the isospin projection T z . Results of fully self-consistent RHB calculations are compared with 
experimental data (for description see the text). 

FIG. 10. Comparison of calculated matter radii for mirror nuclei as function of the isospin 
projection T z . 
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